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The effect of electron-electron correlation on the attoclock experiment 
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We investigate multi-electron effects in strong-field ionization of Helium using a semi-classical 
model that, unlike other commonly used theoretical approaches, takes into account electron-electron 
correlation. Our approach has an additional advantage of allowing to selectively switch off different 
contributions from the parent ion (such as the remaining electron or the nuclear charge) and thereby 
investigate in detail how the final electron angle in the attoclock experiment is influenced by these 
contributions. We find that the bound electron exerts a significant effect on the final electron 
momenta distribution that can, however, be accounted for by an appropriately selected mean field. 
Our results show excellent agreement with other widely used theoretical models done within a single 
active electron approximation. 


The idea of the attoclock with a clearly defined “time- 
zero” for strong field ionization was first introduced in 
[1], and relies on elliptical polarization of an ultra-short 
laser pulse to map the final offset angle of the electron 
momenta distribution to time. Since the attoclock tech¬ 
nique extracts the center of the electron momenta distri¬ 
bution (Fig. 1), the accuracy is limited in principle only 
by the electron statistics at the detector, hence enabling 
the reconstruction of electron dynamics in the attosecond 
(10“^^ sec) domain. However, as many reconstruction 
and experimental calibration procedures in strong field 
ionization experiments mil, the measured delay time re¬ 
lies on a semiclassical (two-step) model within the single 
aetive eleetron approximation IMH]. 

Recent experiments mn\ found an additional angu¬ 
lar offset, relative to an adiabatic two-step model, which 
takes account of the Coulomb field of the parent ion 
within a single active electron approximation (for details 
on the model see 0191EO]). This additional offset was 
explained as due to tunneling delays, following the atto¬ 
clock concept, as presented in [Diziin]. A subsequent 
work [8] showed that using an imaginary time method 
developed by Perelomov, Popov, and Terent’ev (PPT) 
m, which takes account of non-adiabatic effects but 
neglects the Coulomb tail inside the potential barrier, 
gives essentially the same interpretation of the experi¬ 
mental result in [7] as the above mentioned adiabatic 
model. On the numerical front, there remains a dis¬ 
crepancy between time-dependent Schrodinger Equation 
(TDSE) simulations performed within the single active 
electron approximation and the experimental re¬ 

sults in 013 . It has therefore been suggested that this 
discrepancy potentially can be explained by the presence 
of eleetron-eleetron eorrelations in Helium. 

Another recent experiment has shown that electron- 
electron correlations can play an important role in ex¬ 
cited Helium atoms when ionized with an ultra-violet 
photon m- Most recently, however, excellent agree¬ 
ment was obtained for attosecond photoemission mea¬ 
surements from ground state helium with three different 
theoretical models, which go from ab-initio theory tak¬ 


ing into account both electrons to single-active electron 
approximation with both the TDSE and a fully analyt¬ 
ical model m- A question remains as to whether such 
correlation effects are significant in the present experi¬ 
mental range of the attoclock, which involves strong field 
tunnel ionization (using a laser pulse with central wave¬ 
length in the infrared range) of Helium from a ground 
state. Resolving this question with TDSE simulations of 
a two electron atom in the relevant experimental regime 
is a formidable task. We therefore use a semiclassical 
two electron model, introduced in El, that allows for a 
consistent description of the process, taking into account 
the initial bound state and the subsequent propagation 
in the continuum, but, like other semi-classical models, 
neglecting a possible tunneling delay. 

To study the effect of electron-electron interaction on 
the 2-D momentum distribution of the singly ionizing 
electron in He in an elliptically polarized laser field (i.e. 
the experimental observable in 0131111 ), we employ a 
three-dimensional (3-d) semi-classical model for a three- 
body Coulomb system (two electrons plus the nucleus). 
All interactions are fully accounted for including electron- 
electron correlation. The quasiclassical model we use 
entails one electron (electron 1) instantaneously tunnel¬ 
ing through the field-lowered Coulomb potential with a 
quantum tunneling rate given by the Ammosov-Delone- 
Krainov (ADK) formula [H [191 HO]- The longitudinal 
momentum is zero while the transverse one is given by 
a Gaussian distribution m- We compute the exact exit 
point for the tunneling electron using parabolic coordi¬ 
nates m- The remaining electron (electron 2) is mod¬ 
eled by a microcanonical distribution [21] . For the evolu¬ 
tion of the classical trajectories we use the full three-body 
Hamiltonian in the laser field, that is. 
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with E(t) the laser field. An important difference be¬ 
tween our 3-d model and other 3-d semiclassical models 
[TOl [ 22 ] is that our model treats exactly the Coulomb sin¬ 
gularity by introducing regularized coordinates [23| for 
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the time propagation. Transforming to regularized co¬ 
ordinates removes the Coulomb singularity due to the 
electron-nucleus interaction. This results in a faster and 
more stable numerical integration. The model we cur¬ 
rently employ was first described in m- 

Our goal is to compute the offset angle of the 2-D elec¬ 
tron momentum distribution in the plane of polarization 
as a function of the intensity of the laser field and com¬ 
pare with the experimental results in EE]. To do so, we 
model the electric field using a Gaussian pulse, similar to 
the one employed in the experiments that are of interest 
for the current work [SHZlIIHlEll: 

E(t) = Eof(t)(cosc(;tz + esinci;tx) (2) 

with 

f(t) = exp I -21n2 ^i ) (3) 

\ V^fwhm/ J 

and tpwHM = 6 fs, cj = 0.05695 a.u., with Eq the ampli¬ 
tude of the laser field. We use the ellipticity parameter 
e = 0.87. These laser pulse parameters were selected to 
model the experimental set-up described in [6l [7] . The 
use of the dipole approximation in Eq. allows one to 
neglect the spatial dependence of the laser field as well as 
the magnetic fields, which would otherwise considerably 
complicate the dynamics m, hence allowing for an easier 
and more intuitive interpretation of classical trajectories 
[26] . The initial time of propagation is sampled randomly 
from 00 = = [—34,34] a.u., with 0o the initial phase 

of the laser field at the time electron 1 tunnels. 

We compute the 2-D, x-z, momenta distributions of the 
ionizing (tunneling) electron, in single ionization events, 
for different laser field strengths. Specifically, we con¬ 
sider the field strengths 0.03, 0.05, 0.07, 0.09, 0.11 and 
0.13 a.u. In Eig. we plot the 2-D momentum distribu¬ 
tion for a field strength of Eq = 0.09 a.u. We extract the 
offset angle from these momenta distributions as follows. 
Eor each singly ionizing event we register the final mo¬ 
mentum of the tunneling electron, which escapes to the 
continuum. Erom the final momentum we compute the 
angle 0 = arctan(px/pz) and the magnitude of the mo¬ 
mentum p = \/Px + Pz • The distribution S(6>,p) denotes 
the probability for the singly ionizing electron to escape 
with angle 0 and momentum p, respectively. Integrating 
over the magnitude of the momentum, we then obtain 
the distribution f(6>) that denotes the probability for an 
electron to singly ionize with angle 0: 

f(6») = y s(0,p)pdp. (4) 

We use a double gaussian function to fit f(^), i.e. the 
same function as the one used to fit the experimental 
results in ElE]. It was shown in prior work m that, for 
high ellipticities, the radial integration method described 
above integrates exclusively over the transverse momenta 
spreads of the electron distribution at the detector. 


Using the above described procedure, we extract the 
maximum, 6>max, of the f{0) distribution, which corre¬ 
sponds to the most probable value of the angle 0. We 
then identify 90° — ^max as the streaking offset angle (see 
Fig.|^. 90° is how much the momentum at the end of the 
pulse is offset compared to the laser field at the time of 
tunneling, assuming no tunneling delays and neglecting 
the Coulomb interaction with the parent ion. 

Eirst, we investigate whether the electron-electron in¬ 
teraction in the three-body Hamiltonian can be ac¬ 
counted for by using a two-body Coulomb system (one 
electron plus the nucleus) with an effective nuclear charge 
equal to Zeff = 1. In Eig.[^ we plot the 2-D momenta dis¬ 
tributions for Eq = 0.09 a.u. for the three-body and the 
two-body Coulomb systems. The two distributions are 
very similar. Indeed, in Table I we list the offset angles 
90° — 6>max foi* different strengths of the laser field for 
the three-body and the two-body systems. We find that 
the two sets of values are very similar. This similarity 
implies that the electron-electron interaction can be ac¬ 
counted for by a one electron plus a nucleus system with 
an effective charge. 

Next, we investigate further the effect of electron- 
electron correlation on the dynamics of the ionized elec¬ 
tron. Specifically, we compute the offset angle when we 
switch-off the electron-nucleus interaction following tun¬ 
neling. In this case, the offset angle we compute will be 
due to the effect of the Coulomb repulsion that the tun¬ 
neling electron experiences from electron 2 that remains 
bound. In Eig. we show for a fixed field strength that 
as a result of the electron-electron interaction the offset 
angle is in the Pz < 0 & Px < 0 quadrant unlike the offset 
angles for the cases of the full three-body system and the 
two-body system that are in the Pz > 0 & Px < 0 quad¬ 
rant. Specifically, as shown over a wide range of electric 
field strengths in Table II and Eig. while for the full 
three-body and the two-body system the offset angles 
90° — ^max are positive, when we switch-off the electron- 
nucleus interaction on the tunneling electron the offset 
angle is negative (see Eig. 3, dashed red line). This 
implies that the Coulomb repulsion between the bound 
electron and the tunneling electron causes the latter to 
escape to the continuum faster than it would in free-space 
and even more so if we only account for the effect of the 
nucleus on the tunneling electron with Z = 2, (see Table 
II and dashed blue curve in Eig. [^. The above implies 
that the electron-electron interaction (in the mean field 
sense) does have a significant effect on the dynamics of 
the tunneling electron while escaping to the continuum. 
The repulsion from the bound electron causes the tunnel¬ 
ing electron to escape faster compared to the case where 
no electron-electron interaction is present. 

Fig-i shows that the offset angle obtained using a 
semiclassical model (solid red line) is significantly smaller 
than the experimental one, particularly for small inten¬ 
sities. Can electron-electron correlation account for this 
discrepancy? We believe this not to be the case for the 
following reasons: so far we have shown that, within the 
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framework of our model, the offset angle is the same 
whether we fully account for two electron effects or in 
a mean-field sense. This can be seen in Fig.|^ where the 
two electron (three-body system) results (solid red line) 
almost coincide with the single active electron approxi¬ 
mation (solid blue line), when the effective charge of the 
remaining parent ion is such that Zeff = 1. (In regard 
to the single active electron model, the effective charge 
before the ionisation takes place has to be greater than 
one, since the nuclear charge in a ground state of Helium 
is only partially shielded by the other electron.) In ad¬ 
dition, as our results show, electron-electron correlation 
causes the tunneling electron to escape faster to the con¬ 
tinuum resulting in a decreased offset angle. However, 
we note that our theoretical results are obtained with a 
model which does not consider the dynamics inside the 
barrier. It is only after electron 1 exits the barrier that 
we propagate in time accounting for all the interactions 
among the two electrons and the nucleus as well as the 
interactions with the laser field. Fig. [^compares our re¬ 
sults (solid red and blue lines) with experimental data 
and with another adiabatic model described in [22] and 
referred to as TIPIS in [5] (solid black line in Fig. [3). As 
can be seen, there is good agreement between the semi- 
classical models over a wide intensity range in which the 
attoclock experiments in laiT] were performed. 

Investigating electron-electron correlations during the 
tunneling process itself (while the electron is inside the 
barrier) would be the next step in ruling out multi¬ 
electron effects in the attoclock experiments. However, 
such studies are beyond the scope of existing semiclassi- 
cal models. Fully quantum, two-electron, TDSE simula¬ 
tions using a t-SURFF method [28] are currently in the 
planning stages [29]. Such simulations could definitively 
resolve the importance of electron-electron correlations 
in Helium, during strong field ionisation from a ground 
state, by comparing TDSE simulations within a single ac¬ 
tive electron approximation m, with a full two-electron 
TDSE solution. 

In conclusion, we note that taking account of electron- 
electron correlations within a semi-classical model cannot 
account for the additional angular offset measured in the 
recent attoclock experiments Ii[7]. In particular, a care¬ 
ful use of the single active electron approximation results 
in essentially the same angle for the center of the electron 
momenta distribution as the use of a full two electron 
model. One of the advantages of our two-electron model 
is that it allows for a consistent description of the entire 
process, both the initial bound state and the subsequent 
propagation of the ionized electron in the continuum. On 
the other hand, the single active electron approximation 
requires that we use different effective charge for the ini¬ 
tial bound atomic state (where it must be such that it 
reproduces the ionisation potential of Helium) and for 
the propagation in the continuum, where the total effec¬ 
tive charge of the parent ion should be set to one. Using 
the two-electron model also allowed us to investigate the 
relative contributions of the nucleus and the remaining 


Eq = O.OQa.n. 
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FIG. 1. The offset angle for the three-body Coulomb system 
(left) and for the two-body system with an effective charge 
Z = 1 (right). 


bound electron to the final offset angle. 


Eo (a.u.) 

Three-body (°) 

Two-body, Zeft = 1 (°) 

0.03 

8.4 

8.4 

0.05 

6.1 

6.2 

0.07 

5.8 

6.0 

0.09 

6.1 

6.3 

0.11 

6.6 

6.9 

0.13 

6.8 

7.2 


TABLE I. The offset angle for the three-body Coulomb system 
(left) and the two body Coulomb system with an effective 
charge Zeff = 1 (right). 


Eo (a.u.) 

Three-body (a) (°) 

Three-body (b) (°) 

Two-body, Z = 2 (°) 

0.03 

8.4 

-7.1 

17.1 

0.05 

6.1 

-5.5 

12.9 

0.07 

5.8 

-5.2 

11.8 

0.09 

6.1 

-5.2 

12.3 

0.11 

6.6 

-4.9 

12.7 

0.13 

6.8 

-4.6 

13.9 


TABLE 11. The offset angle for the three-body Coulomb prob¬ 
lem as in Table 1 (a), for the three-body problem but with the 
interaction of the tunneling electron with the nucleus switched 
off (b) and for the two-body problem with Z = 2. 
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Eq = 0.09a.u. 
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FIG. 2. The offset angle for the three-body Coulomb system 
when the electron-nucleus interaction is switched off for the 
tunneling electron (left) and for the two-body system with 
Z = 2 (right). 



FIG. 3. Comparison of semiclassical models to experimental 
results. All semiclassical models are within a single active 
electron approximation, except for the solid and dashed red 
lines. The colored dashed lines are the result of selectively 
switching off either the nuclear charge (dashed red line) or 
the charge of the remaining electron (dashed blue line). 
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